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Abstract 

Searches for Dark Matter (DM) particles with indirect detection techniques 
have reached important milestones with the precise measurements of the anti- 
proton (p) and gamma-ray spectra, notably by the Pamela and Fermi-LAT 
experiments. While the 7-ray results have been used to test the thermal Dark 
Matter hypothesis and constrain the Dark Matter annihilation cross section into 
Standard Model (SM) particles, the anti-proton flux measured by the Pamela 
experiment remains relatively unexploited. Here we show that the latter can 
be used to set a constraint on the neutralino-chargino mass difference. To 
illustrate our point we use a Supersymmetric model in which the gauginos are 
light, the sfermions are heavy and the Lightest Supersymmetric Particle (LSP) 
is the neutralino. In this framework the W + W~ production is expected to be 
significant, thus leading to large p and 7-ray fluxes. After determining a generic 
limit on the Dark Matter pair annihilation cross section into W + W~ from the 
p data only, we show that one can constrain scenarios in which the neutralino- 
chargino mass difference is as large as ~ 20 GeV for a mixed neutralino (and 
intermediate choices of the p propagation scheme). This result is consistent 
with the limit obtained by using the Fermi-LAT data. As a result, we can 
safely rule out the pure wino neutralino hypothesis if it is lighter than 450 GeV 
and constitutes all the Dark Matter. 
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1 Introduction 



Indirect searches for Dark Matter, i.e. searches for 'anomalous' features in cosmic rays (e.g. 
gamma-rays, neutrinos, positrons and anti-protons), have been proposed in the late 70's as 
a powerful way to reveal the existence of Dark Matter annihilations in the Milky Way halo 
and beyond [1, 2]. These techniques are meant to give precious insights about the nature 
of the Dark Matter particle and its properties, assuming that a signal is seen. Yet there are 
several limiting factors which weaken their ability to elucidate the dark matter problem. In 
particular indirect detection requires a detailed knowledge of the astrophysical backgrounds 
and foregrounds and therefore depends on the present knowledge of astrophysical sources. 
To make a discovery one either has to carefully remove known (or modelled) background 
in order to expose the 'anomalous' component or hope that the Dark Matter signal is well 
above the background and exhibits very clear features, which would be difficult to mimic 
by invoking astrophysical sources only. 

Currently there are contradicting claims regarding whether indirect detection is giving 
clues of dark matter or not. On one hand, there are possible anomaly detections which could 
be explained in terms of Dark Matter annihilation or decay. These include for example the 
positron excess, as seen in Pamela (and Fermi-LAT) data [3], a possible feature in the 
e + + e~ spectrum [4] , a claimed 7-ray excess at ~ 10 GeV energies [5] 1 and, most recently, 
two possible 7-ray lines at 111 and 129 GeV [13, 14]. On the other hand, a large bulk 
of present astrophysical data essentially seem to validate the modelling of astrophysical 
background sources in the GeV-TeV range (disregarding these possible anomalies), and 
therefore enables one to set powerful constraints on the Dark Matter properties. 

By measuring the gamma-ray spectrum over a large energy range relevant for Dark 
Matter physics, the Fermi-LAT collaboration has been able to set stringent limits on the 
Dark Matter pair annihilation cross section into Standard Model particles. For example, 
using the diffuse 7-ray emission in dwarf spheroidal (dSph) galaxies [15] and also in the 
Milky Way [16, 17], the Fermi-LAT collaboration has ruled out Dark Matter candidates 
with a total annihilation cross section of (av ) = 3 x 10 -26 cm 3 /s if moM J$ 30 GeV. This 
constituted a remarkable milestone as such a value corresponds to that suggested by the 
thermal freeze-out scenario, which is generally considered as a strong argument in favour 
of Weakly Interacting Massive Particles (WIMPs). 

These limits nevertheless weaken at higher DM masses, therefore allowing for heavier 
DM candidates with a larger pair annihilation cross section. For example, for ttt-dm = 
100 GeV the limit relaxes to (av) < 10~ 25 cm 3 /s while for 771dm = 500 GeV, it reads 
(av) < 3 x 10~ 25 cm 3 /s, which is one order of magnitude higher than the 'thermal' cross 
section. 

DM models with such large values of the pair annihilation cross section have actually 
been proposed over the last five years as a consequence of the excesses in e + and e + + 
e~ fluxes. While they may remain hypothetical, discovering such a configuration would 
invalidate the WIMP 'vanilla' model and either point towards the existence of non-thermal 
process in the Early Universe (possibly opening up an unexpected window on fundamental 
physics at high energies) or potentially call for more sophisticated mechanisms, such as 
Freeze-In and regeneration as proposed in [18, 19]. Explaining the observed dark matter 

1 All these claims have possible drawbacks, cf [6, 7, 8, 9, 10, 11, 12]. 
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relic density may remain nevertheless challenging. For example, in [20], it was shown 
that candidates with a total annihilation cross section exceeding (crv) = 10~ 24 cm 3 /s 
(corresponding to a thermal relic density smaller than 3%) would be ruled out by the 
Fermi-LAT experiment if they were regenerated at 100%. 

In addition to measurements of the e + and e + + e~ spectra mentioned above, there is 
also the measurement of the galactic p flux, presented by the Pamela collaboration [21, 
22]. While extensive work was done to explain the electron/positron excesses in terms 
of Dark Matter annihilations (or decays), the implications of the absence of anomalies 
in the p spectrum has remained relatively unexploited. Indeed only a relatively small 
number of works [23, 24, 25, 26, 27, 28, 29] have dealt with it and shown that large Dark 
Matter annihilation cross sections can be constrained by the Pamela data. Among the 
most interesting conclusions which have been reached let us cite for example that in [24] 
constraints on the annihilation cross section into bb were given (for the same mass range 
as is considered in this paper) and limits on the W + W~ final state were mentioned for 
m DM = 1 TeV and one specific set of propagation parameters. In [28], constraints on the 
qqg were set for bino-like neutralinos. 

The first aim of this paper is therefore to propose a more systematic analysis of these 
general anti-proton constraints on the DM annihilation cross section, including paying 
attention to the uncertainties associated with DM and astrophysical predictions. The 
second aim of the present analysis is to demonstrate that these measurements can actually 
constrain the properties of specific DM scenarios, including the mass spectrum in the dark 
sector. To illustrate this, we will work within a 'simplified' version of the phenomenological 
Minimal Supersymmetric Standard Model (pMSSM) [30] in which all sfermion masses are 
set to 2 TeV, except for the stop and sbottom masses. The soft masses for the stop 
are allowed to be much lighter to obtain a Higgs at 125 GeV. In this scenario the only 
particles with masses below the TeV threshold are therefore the neutralino, chargino, the 
supersymmetric Higgses and the lightest stop and sbottom. Such a configuration of 'light' 
gauginos and heavy sfermions may actually seem unnatural from a supersymmetric point 
of view (albeit close to split SUSY [31]) but it is supported by the unfruitful searches for 
squarks and gluinos at LHC, at least to some extent 2 . 

With this very set up in mind, one can investigate scenarios where the neutralino pair 
annihilation cross section into W + W~ gauge bosons is enhanced (due in particular to the 
chargino exchange diagram). Such a large annihilation cross section gives both a significant 
anti-proton and diffuse gamma ray flux, together with a gamma ray line, and is therefore 
potentially constrained by the Pamela and Fermi-LAT data. In Supersymmetry, such an 
enhancement is realised when the LSP neutralino is mass degenerated with the chargino, 
i.e. when the neutralino has a significant wino component. The combination of both 
Fermi-LAT and Pamela data is therefore expected to constrain the wino fraction of the 
lightest neutralino, thus realizing our second aim. Note that constraints on the neutralino 
composition are also expected to be obtained in presence of a lower sfermion mass spectrum. 
However the effect of the chargino-neutralino mass degeneracy on 7— ray and p production 
would be much harder to characterise. Hence our choice in favour of a heavy sfermion mass 
spectrum. 

2 Even though, admittedly, those negative searches may also be a sign that Supersymmetry is not 
realised at the TeV scale. 
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Table 1: Propagation parameters for anti-protons in the galactic halo (from [32, 33]). Here 
S and /Co are the index and the normalization of the diffusion coefficient, V conv is the velocity of 
the convective wind and L is the thickness of the diffusive cylinder. 

The paper is organised as follows. In section 2 we derive generic constraints on the 
Dark Matter pair annihilation cross section into W + W~ from anti-proton data and recall 
the Fermi-LAT limits that are obtained from gamma-ray observations in the Milky Way 
and dwarf Spheroidal galaxies. In Section 3 we present the Supersymmetric model that we 
shall consider and explain how we perform the scans of the parameter space. Finally in 
Section 4 we apply the Pamela and Fermi-LAT limits to our SUSY model and show that 
the anti-proton data can be more constraining than gamma-ray observations. We conclude 
in Section 5. 

2 Anti-proton and 7-ray bounds on <t D m dm w+w~ 

In this section we discuss how anti-proton and gamma ray data impose generic constraints 
on the Dark Matter pair annihilation cross section into W + W~ as a function of the Dark 
Matter mass. 

2.1 Generic bounds on ct D m dm w+w- from anti-protons 

W ± production in space leads to abundant anti-proton production as the W ±, s decay 
products hadronize. The flux of anti-protons thus produced by DM annihilations into a 
pair of W ± gauge bosons in the Milky Way and collected at Earth is therefore determined 
by the Dark Matter pair annihilation cross section into W + W~, the Dark Matter mass and 
the Dark Matter halo profile. It also depends on the anti-proton propagation parameters 
which are being considered. Hereafter we will assume that the dark matter halo profile is 
well described by an Einasto profile (we checked that other choices make a small difference) 
and consider the standard three sets of propagation parameters ('MIN', 'MED', 'MAX') 
summarised in table 1. In practice, we use the anti-protons fluxes which are given in [34], 
to which we refer for further details. 

In order to constrain the annihilation cross section, we will consider that all present data 
define the maximal flux in anti-proton that is allowed by the PAMELA [22] experiment 3 . 
Both the predicted energy spectrum and the flux depend on the dark matter mass that is 
being assumed. For each value wdm, we will therefore compare the sum of the astrophysical 
background flux and predicted anti-protons flux originating from Dark Matter with the 
Pamela data. Given the uncertainties on the astrophysical background, we will apply 
two different procedures to derive meaningful limits. One can be regarded as aggressive 

3 To avoid the uncertainty related to solar modulation, we restrict ourselves to using the Pamela data 
above an anti-proton energy of 10 GeV. 
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(it assumes a fixed background) while the other one is more conservative (the background 
can be adjusted within the uncertainties). 

o For obtaining aggressive limits (referred to as fixed background in the following), we 
adopt the standard flux of astrophysical (secondary) anti-protons from [35] and add 
it to the DM anti-protons flux. We then compare the result with the Pamela data 
and derive a 95% CL. limit by imposing that the global x 2 of the background + DM 
flux does not exceed by more than 4 units the \ 2 of the null hypothesis (background 
only). 

o For obtaining conservative limits (hereafter referred to as marginalized background) , 
we take again the standard form of the background spectrum predicted in [35], except 
that now we allow for the normalisation of the background spectrum A and the 
spectral index p to vary within 40% and ±0.1 respectively (for each value of the DM 
mass and pair annihilation cross section into W + W~). 

In practice, we multiply the standard description of the background spectrum by a 
factor A (T/T ) p , where T is the anti-proton kinetic energy, T = 30 GeV is a pivot 
energy and with 0.6 < A < 1.4 and —0.1 < p < +0.1. These are quite generous 
intervals, which allow to include the uncertainty predicted in [35]. We then add 
up the DM contribution expected for each point in the parameter space defined as 
(rriDMi ( av )) an d identify the pair of parameters A and p which minimises the global 
X 2 with the Pamela data. This procedure therefore corresponds to marginalising 
over the parameters of the uncertain astrophysical background point-by-point in the 
DM parameter space. Again, the 95% CL. is then imposed by requiring that the 
marginalised global \ 2 does not exceed 4 units with respect to the null hypothesis 
(which has been marginalised consistently). 

By considering a variable background spectrum (within the uncertainties) for each 
value of the DM mass and cross section, we can increase the gap between the expected 
p background and the actual Pamela data. As a result this leaves more space for 
a possible DM injection of anti-protons and leads to weaker limits on the DM pair 
annihilation cross section. A similar approach was used in [36] but to reduce the gap 
between the astrophysical background and the data. 

Our constraints are displayed in Fig. 1 for the 'MIN','MED', 'MAX' set of parameters. 
As expected, the 'conservative' limits are slightly less constraining than the 'aggressive' 
ones. Also we find that the choice of propagation parameters has a big impact on the type 
of constraints that can be set: in terms of cross sections, the difference between the 'MIN' 
and 'MAX' limits exceeds a factor 10. 

To understand more precisely how these constraints work, we defined 5 scenarios (here- 
after referred to as 'A','B','C','D','E'), corresponding to different DM masses, cross sections, 
propagation parameters and constraint procedures. The corresponding fluxes are plotted 
in Fig. 2. As one can see, benchmark points 'A' to 'D' correspond to 'borderline' scenar- 
ios where the total p flux (i.e. the sum of the expected flux from DM and astrophysical 
background) is not significantly exceeding the data. Point 'E', on the other hand, displays 
'how badly' the data is violated inside the excluded region. 

The first apparent feature from Fig. 2 is that one can actually exclude a small excess 
in anti-protons produced by relatively light Dark Matter particles because the Pamela 
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Figure 1: Anti-proton constraints on 

DM annihilation into W + W~ . The up- 
per left, upper right and lower panels 
refer respectively to the 'MIN', 'MED' 
and 'MAX' propagation parameters. The 
constraints obtained by the Fermi-Z^4T 
collaboration from satellite dwarf galax- 
ies are superimposed. We also display 
five benchmark points. 



data set have very small error bars at energies below 100 GeV, hence the strength of the 
constraints. It is then instructive to compare case 'A' and 'B': these two scenarios refer to 
the same DM mass and constraint procedure; they also predict a very similar flux, as can be 
seen in Fig. 2, but have a different annihilation cross section. The latter is much larger for 
'A' than for 'B'. This is because the propagation scheme was assumed to be 'MIN' for the 
former and 'MED' for the latter. With the 'MIN' propagation set, the yield of anti-protons 
is about one order of magnitude smaller than with 'MED' (since the galactic diffusion zone 
is much smaller in the former case) and therefore the constraint on the annihilation cross 
section is about one order of magnitude looser than for the 'MED' case. On the other hand 
the constraint obtained for 'MAX' (which is not shown here) is stronger than for 'MED'. 

The comparison between points 'B' and 'C shows the impact of the constraint proce- 
dures. Although both 'B' and 'C have the same DM mass and propagation scheme, we 
find that the value of the annihilation cross section that is allowed for 'C is larger than 
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for 'B'. The reason is that ! C corresponds to the scenario in which the limit is obtained 
by using the 'marginalised background' procedure (i.e. where the background is allowed 
to retract within the uncertainties) so there is more room for DM while 'B' corresponds to 
a 'fixed background' scenario so the associated constraints are stronger. 

Finally, the comparison between 'C and 'D' enables one to understand why the 'margi- 
nalised background' constraints are rather independent of the DM mass, despite the fact 
that the error bars in the Pamela data become larger at larger energies. For a large DM 
mass (case 'D') the p flux is shifted towards larger energies and rather negligible at ~ 10 
GeV with respect to the astrophysical background; there is thus little room to reduce the 
the bakcground (which alone has to fit the data at low energy) and consequently there is 
little room left for a DM contribution at large energies. As a result, the bound remains 
stringent. 

2.2 Generic bounds on ct D m dm^w+w- from gamma-rays 

In DM scenarios, the W production is associated with gamma-ray emission through (i) 
the decay and hadronisation of the decay products of the W ± bosons, (ii) the radiation 
of a photon from the internal and/or final states associated with DM DM — > W + W~ (iii) 
DM annihilations into 77 and 7Z (which can be seen as a higher order process based on 
DM DM — > W + W~). The first case leads to a continuum spectrum of 7-rays (the energy 
spectra can be e.g. found in [34], for any value of the DM mass); the second leads to sharp 
features in the 7-ray continuum spectrum and the third to 7-ray lines. The resulting fluxes 
from these process have to be compared with the gamma-ray flux measurements from the 



7 



Milky Way or from other nearby galaxies. Therefore we now review the current 7-ray 
constraints derived in the literature (mainly from Fermi-LAT analyses), paying particular 
attention to that derived from the W + W~ channel since this is the main focus of our 
analysis. 

2.2.1 Continuum 

The Fermi-LAT collaboration has recently published two different analyses of the contin- 
uum diffuse gamma-ray emission from the Milky Way halo [16, 17]. Since no clear DM 
signal has been found, these have been used to set upper limits on the DM pair annihilation 
cross-section into various channels: e.g. bb, gg, W + W~, e + e", /i + /U~, t + t~. For relatively 
light DM (~ 20 GeV) and e.g. the bb channel the limits reach the canonical value of the 
cross section (namely (av) = 3 x 10~ 26 cm 3 /s), provided that the most aggressive proce- 
dure is used. For DM masses (9(100) GeV and for the W + W~ channel the limit reads 
(av) <2x 10~ 24 cm 3 /s. However the most stringent limits on the DM annihilation cross 
section have actually been obtained from another Fermi-LAT analysis based on the diffuse 
7-ray emission from dSph galaxies; these Dark Matter dominated objects indeed represent 
a good target for Dark Matter searches. 

In the present analysis we will use the results from [15, 37] (see also [38]). Although 
they use slightly different sets of targets 4 , slightly different datasets 5 and a different 
analysis procedure ([37] introduces a frequentist Neyman construction), they both derive 
consistent limits for the bb, W + W~ , fi + ,t + t~ channels. If we apply - for definiteness 
- the constraints from [15] and assume a DM mass of 100 GeV, the limit for the W + W~ 
channel reads as (av) < 8.5 x 10~ 26 cm 3 /s. The analysis procedure in [15] allows one to 
incorporate the uncertainties associated with the DM energy density profile of individual 
dSph galaxies, which was shown to lead to an error band of about an order of magnitude 
on the constraint in [37]. Here we do not attempt to address these issues; we simply draw 
the attention of the reader that these constraints have to be taken with care until a better 
determination of the DM energy density profile in dSph galaxies is available. Consequently, 
we adopt the rather conservative constraints in this paper. 

In Fig. 1 we compare the dSph galaxies limits with the Pamela anti-proton bounds 
that were derived in Sec. 2.1. We see that, depending on the propagation scheme that 
has been chosen for the anti-protons, the dSph galaxies 7-ray bounds is somewhat more 
stringent or looser than the constraints from the anti-proton data. For example, for the 
'MED' case and 'marginalized background', the p limits becomes more constraining than 
the 7-ray bounds when todm > 290 GeV. However they are stronger than the 7-ray limits 
whatever the value of moM (assuming mDM > 100 GeV) for a 'fixed' background. Since 
nevertheless the p and 7-ray limits are basically of the same order of magnitude, we will 
include both constraints in our study. 

4 Ref. [37] uses 7 dSphs —Bootes I, Draco, Fornax, Sculptor, Sextans, Ursa Minor, and Segue 1 while 
Ref. [15] uses 10 dSphs —the same as above plus Carina, Coma Berenices and Ursa Major II—. 

5 Ref. [15] uses 24 months between August 2009 and August 2010 while Ref [37]) uses 3 years between 
August 2009 and August 2011. 
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2.2.2 Internal Bremsstrahlung and Final State Radiation 

Gamma rays produced directly as radiation from an internal line or a final state are in 
general suppressed by the fine structure constant, a. However, for a t-channel diagram, 
the associated cross section can be enhanced when the intermediate particle is almost 
mass degenerated with the DM. Typically the enhancement factor is about mp M / (M| — 
m Divi) where Mj the mass of the intermediate particle (i.e. a chargino for neutralino 
pair annihilation into a W pair). These process are model dependent and cannot be 
constrained generically but they will be included in our 7-ray estimates when we investigate 
the neutralino pair annihilations into W + W~ in the pMSSM. 



2.2.3 Line(s) 



Annihilations directly into 77 or 7Z occur at one-loop level (since DM particles do not 
couple directly to photons) and are therefore generically suppressed. However they lead to 
a distinctive signature, namely a mono-energetic gamma-ray line at an energy E = jbdm 
or E = wdm (1 — m z/(^ m DM)) which can be looked for. 

With possible evidence for two gamma-ray lines at 129 and 111 GeV (which have been 
speculated as originating from DM particles with a mass of about 130 GeV annihilating 
into 77 and jZ), indirect detection of DM particles seem promising. Yet the existence 
of these lines remain to be confirmed by the Fermi-LAT collaboration and their origin 
to be shown as being exotic. Since the purpose of this study is to set constraints on the 
DM properties (and owing to these uncertainties on the existence and origin of these lines) 
we will disregard the results of [13] and only consider the constraints which were reported 
by the Fermi-LAT collaboration on line searches in the Milky Way [17], where the upper 
limits on cxi;dmdm->77 and 0"Wdmdm-s>Z7 range from 0.03 to 4.6 x 10 _27 cm 3 s _1 and 1 to 
10 x 10~ 27 cm 3 s _1 respectively, for DM masses up to 200 GeV. Constraints on <7fDMDM-m 
were also obtained from dSph galaxies [39] but they are not as stringent as those obtained 
from the Milky Way. 

Since the status of these searches is not definite, we made the choice to not include 
these constraints to perform the scans over the pMSSM parameter space. However we do 
check that the scenarios which survive the p and 7-ray constraints are not killed by these 
line searches. 
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Figure 3: Dominant neutralino pair annihilation diagrams into W + W ,77 and 7 Z for 

this analysis. 
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3 Chargino-neutralino mass degeneracy 



Now that we have obtained the maximal value of the Dark Matter pair annihilation cross 
section into W + W~ that is observationally allowed as a function of the Dark Matter mass, 
we can focus on a specific Supersymmetric model and investigate the impact of this generic 
limit on the neutralino Dark Matter parameter space. 

3.1 Neutralino pair annihilations into W + W~ 

In a scenario where all the sfermions are very heavy, the dominant neutralino annihila- 
tion channels are expected to be mostly into gauge bosons, more specifically into W + W~ 
pairs. All loop-induced W production diagrams which involve sfermions are expected to 
be suppressed. Hence the process which are expected to lead to a significant W + W" pro- 
duction in the pMSSM only involve charginos and Z boson. The corresponding diagrams 
are displayed in Fig. 3. Since they correspond to s— and t— channel diagrams, we typically 
expect resonant or enhanced annihilations when m x o ~ mz/2 or m x o ~ m x ± (i.e. when the 
neutralino and chargino are mass degenerated). These ultimately enhance the neutralino 
pair annihilations into 77 [40, 41] and 7Z [41, 42] through in particular the two 'loop' 
diagrams displayed in Fig. 3. 

The questions that we want to address in the next subsections are: i) which part of 
the SUSY parameter space is excluded by the p limits and do these limits exclude more 
allowed configurations than the 7-rays bounds? ii) which values of the neutralino-chargino 
mass degeneracy are actually constrained by astrophysical data? 

3.2 Exploring the supersymmetric parameter space 

To answer this, we will explore the pMSSM parameter space using the same Markov Chain 
Monte Carlo method as in [43] coupled to the micrOMEGAs code [44] and the SOFTSUSY 
spectrum calculator [45]. 

Our free parameters and their corresponding range are summarised in Table (2). These 
include the soft mass terms associated with the squarks of the third generation (i.e. Mg 3 
and Mu 3 ) and the trilinear coupling A t . To obtain sfermion masses at the TeV scale, we 
set all the soft masses to 2 TeV. In addition, we set the trilinear couplings to and the 
CP-odd Higgs mass to 1 TeV. In this framework, the bino mass Mi does not exceed 500 
GeV; our choice for the other parameters indeed ensures that the neutralinos and charginos 
are light and the mass splitting between the neutralinos and charginos remains relatively 
small. 

On top of these free parameters, we had to include some nuisance parameters over 
which we will marginalise [46]. These are related in particular to the quark content of the 
nucleons (since they have a non-negligible impact on the computation of the Dark Matter- 
nucleon scattering cross section) and the top mass (since it has an impact on the Higgs 
sector). All of them are allowed to vary in the range [N exp -3a, N exp +3cr], with N exp (a) 
the corresponding experimental value (error), as shown in Table 2. 

We also require that the lightest Higgs mass only varies within the range allowed by 
the ATLAS and CMS experiments [47, 48], namely = 125.9 ± 2.0 GeV. However, by 
precaution, we checked that the scenarios which seemed allowed were compatible with the 
latest version of the HiggsBounds code [49] (even though the most recent LHC results on the 
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O'-kN 

m t 


0.553 ± 0.043 [51] 

18.9 ± 0.8 [51] 
44 ± 5 MeV [52] 
21 ± 7 MeV [52] 
173.2 ± 0.9 GeV [53] 


C^mjrrid, 0.51, 0.596, 0.043) 
£i(m s /m d ,18.1,19,7,0.8) 
£1(^,39,49,5) 
£i((T s , 14, 28,7) 
£i(m t , 172.3, 174.1, 0.9) 



Table 2: Range chosen for the pMSSM free parameters and nuisance parameters. 



Higgs [47, 48] are not included in this version). Note that we did not add any requirement 
about the Higgs signal strength to perform the scans. Would ATLAS and CMS confirm an 
'anomalous' Higgs signal strength into 7 7 (i.e. larger than SM expectations) with a high 
confidence level, the pMSSM would be difficult to reconcile with the data. However the 
principles of our analysis would remain valid and could still be used to constrain small mass 
degeneracies between the Dark Matter and another (e.g. t— channel exchange) intermediate 
particle. 

The neutralino relic density is allowed to vary between Q x o h 2 £ [1% WMAP7, WMAP7] 
with n WMAP7 h 2 = 0.1123 ± 0.0035, using WMAP 7-year + BAO + H and the RECFAST 
1.4.2 code [50] . We do not consider smaller values of the relic density as these correspond 
to DM scenarios with very large values of the annihilation cross section and ultimately 
overproduce gamma-rays in the galaxy (i.e. are excluded) if their relic density is entirely 
regenerated, cf [20]. 

For each scenario (corresponding to a point in the pMSSM parameter space), we then 
calculate the total likelihood function. The latter is a product of likelihood functions 
associated with each observable, nuisance parameters and free parameters which have been 
chosen according to the criteria described below. 



3.2.1 Z2i(x, x m j n , x max , <t) 

To m h , Q x oh 2 and all nuisance parameters, we associate a likelihood function £ x which 
decays exponentially at the edges of a well-defined range [x m in, x max \ with a variance a: 

if X <C ^mm? 




C\(yX^ X m i n , X max , Oj \ - ^2 ~ \f T ~> T ("0 



Here x is either m/,, the LSP relic density or the nuisance parameters. Note that we assume 
flat prior for all nuisance parameters. For the free parameters, we will consider a slight 
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modification of the above function, namely 



£/ \ JO for x <C x m i n or x !> x maxj , > 

1\X, X m i n , X max , U ) Sir r i v / 

I 1 tor x £ [x m i n , x max \ 

so as to immediately reject all the scenarios in which one of the free parameters falls outside 
of the specified range. In fact we also immediately reject points where the neutralino is not 
the LSP, where the LEP limits on chargino, slepton and squark masses are not satisfied or 
the calculations of the SUSY spectrum fail. We did not implement LHC limits on sfermion 
masses because our requirement of a sfermion spectrum at the TeV scale should ensure 
that they are satisfied. However updates on direct searches for relatively light stop and 
sbottom would be useful to implement to further constrain the parameter space. 



• >._._ L*2 \X, X m i n , X max , CT) 

We will use a Gaussian Likelihood function, £ 2 , for the B(b — > X*j) observable (one of the 
B-physics observables that we consider) with experimental mean value \i and theoretical 
+ experimental error a : 

_ O-m) 2 

C 2 (x,fi,cr) = e 2« 2 . (3) 

These observables are important as they receive a potentially large contribution from 
chargino/stop loops when either one of these particles is light. This contribution can be 
compensated by the charged Higgs/top diagram but the latter is however suppressed when 
the charged Higgs mass is at the TeV scale. 



O.Ji.O Z_-3^X, X m i n , X max , <T ) 

We also include a Likelihood function £3 for the 2012 XENON100 limits [54] to ensure that 
the scans do not select too large values of the Dark Matter-nucleon scattering cross section. 
In fact we also associate £3(0;, /i, a) to regions of the parameter space where av x o x o_+ w + w - 
is greater than 10~ 27 cm 3 /s. The latter is defined as follows: 

C 3 (x,fi,a) = \^=£- (4) 

1 + e ° 

where the lower or upper experimental bound are associated with the positive or negative 
variance a respectively. Note that some experimental measurements are very discrepant 
with the SM expectations (namely the anomalous magnetic moment of the muon Aa M and 
the branching ratio B(B + — > t + v t )). These observables receive additional contributions 
from particles in the pMSSM but they are too small to explain the observations. Therefore 
we associate a Likelihood function to them which corresponds to Cs(x,fi,a) so that the 
Likelihood is equal to unity if the predictions are much below the measured value. The set 
of constraints that we use is summarised in Table 3. 



4 Results 

The results of our scans are shown in Fig. 4. In the upper left panel is displayed the 
neutralino pair annihilation cross section into W + W~ as a function of the mass degeneracy 
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Constraint 


Value /Range 


Tolerance 


Likelihood 


m h (GeV) [47, 48] 


[123.9, 127.9] 


0.1 


£i(m h , 123.9, 127.9, 0.1) 


n v oh 2 [501 


[0.001123, 0.1123] 


0.0035 


£i(fi x o/i 2 , 0.001123, 0.1123, 0.0035) 


B(b -> X*7) x 10 4 

V Si) 


3.55 


exp : 0.24, 0.09 


£ 2 (10 4 S(6^ X* 7 ),3.55, 


[55, 561 




th : 0.23 


V0.24 2 + 0.09 2 + 0.23 2 ) 


a s L (pb) 


(mnM, C/v) plane 
from [54] 


<7/v(mr)M)/100 


'Csfo" s 'o v , (7/v(r7inM), — c r /v(mnM)/100) 


(TV O.w+w — 
X^X]_ — ^ vv ~ W 


1 


0.01 


C-\{ov o o_ i T^4- , 1, 0.01 1 


(10~ 27 cm 3 /s) 








Aa M x 10 10 [57] 


28.70 


0.287 


£ 3 (10 10 Aa Al , 28.70, -0.287) 


B(B S -> m + aO x 10 9 [58] 


4.5 


0.045 


r 3 (10 9 S(B S -> /x + M"), 4.5, -0.045) 




0.002 


0.0001 


£ 3 (Ap, 0.002, -0.0001) 


#B + ^r+P T ( P SM ) t 59 ] 


2.219 


2.219xl0- 2 


£ 3 (-Rb+-*-+p t , 2.219, -2.219 x 10~ 2 ) 


Z x?X? (MeV) 


1.7 


0.3 


£ 3 (^^X?X?,l-7,-0.3) 




1 


0.01 


^3(o- e +e-^-x?xE,3 X 


S(X§, 3 ^X?) (Pb) [60] 






iB(4 3 ->ZxS), 1,-0.01) 



Table 3: Constraints imposed in the MCMC, from [61] unless noted otherwise. 



between the neutralino and the chargino and in terms of the neutralino composition. In 
the upper and lower right panels we show the pair annihilation cross section into 7Z and 
77 respectively as a function of the neutralino-chargino mass degeneracy Am = m x + — m x o 
and in the lower left panel we give the forecasted spin-independent elastic scattering cross 
section as a function of the neutralino mass for a Xenon-based experiment. 

The left upper panel indicates the neutralino composition which maximises the W ± 
production. As one can see scenarios where av x o x o^ w+w - is the largest and the neutralino- 
chargino mass splitting is the smallest correspond to neutralinos with a very large wino 
fraction. Large values of both av x » x o_^ w + w - and the Xi ~ X + mass splitting correspond on 
the other hand to wino-dominated neutralinos but with a non negligible higgsino compo- 
nent. For these two types of wino-dominated configurations the neutralino and chargino 
mass degeneracy is large enough to make the t— channel (chargino) exchange diagram very 
large. As the wino fraction decreases, the mass splitting becomes larger and the t— channel 
chargino exchange diagram contribution decreases. However it remains large till the hig- 
gsino fraction which ensures large values of the Xi~X + ~ W~ coupling remains significant 
(i.e. dominates over the bino fraction). 

The upper right panel of Fig. 4 shows which values of the neutralino pair annihilation 
cross section into Z7 are excluded by astrophysical data as a function of the neutralino- 
chargino mass splitting. A similar plot is shown for 77 but the colour code now illustrates 
the relation between the different values of this cross section and the neutralino 'thermal' 
relic density. As one can see the shape of the scenario distribution for 77 and is 
essentially the same in the (Am, av) plane. However the Z7 cross section is approximately 
10 times larger than that for 77 for every scenario. Hence combining these two figures 
actually gives an information about the relic density of the scenarios which are excluded 
by astrophysical data. 

In the Z7 plot (upper right panel of Fig. 4), the points excluded by the Fermi-LAT 
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dSph continuum 7-ray data are displayed in yellow (we do not superimpose the constraints 
from line searches). Those correspond, by construction, to scenarios where there is a very 
large W ± production (and thus a large contribution to the continuum 7-ray spectrum) but 
also to the few regions in which the LSP is heavy and where the bb final state (associated 
with the s— channel pseudo-scalar Higgs exchange and which cannot be discarded as it is 
significant) overproduces 7-rays. The regions which are excluded by the Pamela data are 
shown in red. The black points correspond to scenarios excluded by both the Pamela and 
Fermi-LAT data while those in green represent the points allowed by these two types of 
constraints. 

As one can see from the distribution of black points the largest values of the annihilation 
cross sections into Z7 (and therefore jy + jy~) are excluded by both measurements. Since 
these scenarios correspond to a small (or relatively small) chargino-neutralino mass splitting 
and thus large values of the t— channel chargino exchange diagram, we can conclude that 
both Pamela and Fermi-LAT data are relevant to constrain wino-dominated neutralinos. 
A small number of these configurations is however constrained by only one of the Pamela 
or Fermi-LAT dataset but this does not affect the maximal value of the Xi ~~ X + mass 
splitting that can be excluded by using astrophysical considerations. 

By inspecting where the neutralino pair annihilations into Zj, 77 and W + W~ are 
significant in these plots, one also finds that higgsino-dominated scenarios are constrained 
by both Pamela and Fermi-LAT data because the box diagram (cf the lower left diagram 
in Fig. 3) still generates a large W ± production. In fact, for such a LSP, the annihilation 
cross section into ZZ also becomes non-negligible compared to that into W + W~ . Since the 
expected 7-ray and p spectra from W ± and Z production are very similar, we accounted 
for them both when we made the comparison with the PAMELA and Fermi-LAT data. 

Finally the green points which pass all the constraints have a non-negligible bino com- 
ponent. This reduces the chargino exchange diagram contribution and thus enables to 
decrease the W ± (and therefore anti-proton and 7-ray) production. For these bino-like con- 
figurations one expects the stop and sbottom exchange to be relevant, leading to quarks 
in the final state and possibly (in particular for bb) an overproduction of gamma-rays. 
Note that such process would also compete with the neutralino pair annihilation into SM 
fermions near pseudo-scalar Higgs resonances for heavy neutralinos. 

The lower left panel of Fig. 4 indicates whether the spin-independent elastic scattering 
cross section is compatible with the latest results from the XENON100 experiment [54]. 
Again in green are the points which are astrophysically allowed, in black the points which 
are excluded by both Pamela and Fermi-LAT data and in red or yellow the points which 
are either excluded by the Pamela or Fermi-LAT experiments respectively. Clearly one 
can see that the combination of both the Pamela and Fermi-LAT astrophysical con- 
straints surpass the latest exclusion limit set by the XENON100 experiment. In fact in 
general the astrophysical constraints discussed in this paper even have a stronger exclusion 
power than the forecasted XENON1T limit, illustrating how important adding astrophys- 
ical knowledge is in this specific scenario. 

Even though many configurations are excluded by the Pamela and Fermi-LAT data, 
we do find scenarios which are neither excluded by the XENON100 2012 limit nor by the 
astrophysical constraints discussed in this paper. Hence the XENON100 experiment could 
still discover evidence for relatively light pMSSM neutralinos {m x o < 500 GeV) if these 
particles indeed exist. We note nevertheless that in [62], a constraint as strong as the 
XENON100 2012 limit was obtained by using the XENON100 2011 data and a Bayesian 
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Figure 4: Plots of the neutralino pair annihilation cross section into W + W~ (left upper panel) 
and "fZ (right upper panel) as a function of the chargino-neutralino mass splitting and the spin- 
independent DM-nucleon cross section as a function of the Dark Matter mass (lower left panel) 
together with the XENON 2012 limit. The Freeze-Out relic density is displayed in the lower right 
panel for the annihilation cross section into 77 as a function of the mass splitting. 



analysis where the full information available in the (Si, S2) scintillation plane was exploited. 
It is therefore likely that the XENON100 experiment can improve its present exclusion limit 
with the 2012 data and rule out some of the configurations shown here in green. 

In these figures we have assumed that the relic density was regenerated at 100 % for 
candidates with a total annihilation cross section much larger than the 'thermal' one (i.e. 
with a suppressed Freeze-Out relic density). This way we could ensure a fair comparison 
between theoretical expectations and the limits set by the Fermi-LAT and XENON100 ex- 
periments. Looking at the av x o x o_^. 71 plot, one sees that invoking regeneration is needed for 
all scenarios with a chargino-neutralino mass splitting smaller than ~ 20 GeV 6 . Assuming 
that all these candidates have the correct relic density, we could indeed exclude scenarios 
with a neutralino-chargino mass splitting up to 20 GeV and values of av x o x o^ yZ down to 
10 -28 cm 3 /s (see Fig.4), corresponding to av x o x o^ w+w - > 10~ 25 cm 3 /s and Vth 2 <C 0.06. 
However, relaxing the regeneration assumption would completely relax the exclusion re- 
gions and therefore the bound on the mass splitting (apart perhaps from scenarios with 
extremely small mass splitting). 

As a side comment regarding the so-called '130 GeV line': we do find scenarios where 
<7 " l 'xix -+77 ~ 10~ 27 cm 3 /s, which is the value of the cross section that is required to explain 
the feature in the spectrum. These configurations predict a neutralino-chargino mass 

6 For larger values of the mass splitting, no regeneration assumption is required but the annihilation 
cross sections into 77 and 7Z are strongly suppressed. In particular cv X iXi^ii * s muc h below 10 -29 cm /s. 
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Figure 5: Annihilation cross section into W + W as a function of the neutralino Dark 
Matter mass. The chargino NLSF '-neutralino LSP mass splitting is shown as colour code. 

splitting greater than ~ 0.2 GeV. However none of the points corresponding to neutralinos 
with a mass of about 130 GeV are allowed by the Pamela data. Hence, our results 
suggest that one cannot explain the '130 GeV line' in our simplified version of the pMSSM, 
which is in agreement with [63, 64]. Indeed, due to the anti-proton limit, scenarios with 
cxf x o x o^ 77 ~ 10~ 27 cm 3 /s rather correspond to neutralinos with a mass of about 400- 
450 GeV. In fact, for the same reason, all the points with o~v x o x o_+ yj > 2 x 10~ 28 cm 3 /s 
correspond to configurations where m x o >200 GeV. Finally note that in the pMSSM the 
existence of 130 GeV neutralinos should give rise to a second 7-ray line at ~ 111 GeV (on 
top of that at 130 GeV), corresponding to the neutralino pair annihilation into 7Z. Given 
our prediction for 7Z and 77, the flux associated with this 111 GeV line should be about 
ten times larger than that corresponding to the 130 GeV line, which is in conflict with the 
observations. 

Finally, in Fig.5, we show the annihilation cross section into W + W~ as a function 
of the neutralino mass and superimpose the Pamela (for the 'MED' set of propagation 
parameters and marginalised background, i.e. the conservative limits) and Fermi-LAT 
limits (red and yellow lines respectively). The colour code indicates the different values of 
the neutralino-chargino mass splitting. As can be seen from this plot, the Pamela and 
Fermi-LAT constraints are actually complementary. The Fermi-LAT limit excludes more 
configurations below 300 GeV than the Pamela bound but it assumes that the observations 
are independent of the Dark Matter energy distribution in dSph galaxies, which can be 
debated [65, 66]. In contrast, the anti-proton limit excludes a bit more configurations than 
the gamma-rays above 300 GeV. This is reassuring since it is set by observations 'within' 



16 



the galaxy but the drawback is that it relies on a specific choice of propagation parameters 
and astrophysical knowledge of astrophysical sources. In any case, the fact that both limits 
exclude similar configurations enables us to validate the exclusion region that we found. 

Hence the main information that one can read of from this plot, combined with that 
displayed in Fig. 4, is that: 

• one can rule out neutralino-chargino mass splitting up to ~ 20 GeV if m x o < 150 
GeV and the neutralino is a mixture of wino and higgsino 

• one can exclude all scenarios in which the wino-chargino mass difference is smaller 
than 0.2 GeV for m x o < 500 GeV, thanks to both Pamela and Fermi-LAT data. 

5 Conclusions 

In this paper we explicitly derived the constraints on the DMDM — » W + W~ annihilation 
cross section by using the Pamela anti-proton data and paying particular attention to 
the choice of propagation parameters and uncertainties on the astrophysical background. 
Our results are independent of the so-called Pamela positron excess and are obtained for 
two different (fixed vs marginalised) choices of the background spectrum; they are also 
consistent with the enhancement factor which was derived in [24] and the detailed analysis 
of [26], for the cases where the propagation parameters overlap. 

We then compared these bounds with the most stringent gamma-ray limits which have 
been derived using the Fermi-LAT measurements of the gamma-ray continuum spectrum 
from dSph galaxies, for the same DM annihilation channel and DM mass range. We 
found that the anti-proton constraints appear to be very competitive with the gamma- 
ray bounds. More precisely, choosing the 'MED' propagation scheme, the p constraints 
are slightly weaker than the 7-ray ones when moM < 300 GeV and slightly stronger for 
mDM > 300 GeV. On the other hand, the anti-proton constraints are stronger if we assume 
the 'MAX' set of propagation parameters and less powerful if we assume the 'MIN' set. 
We also recall that the gamma ray limits themselves may be subject to some uncertainties 
related to the modelling of the DM profile in dSph galaxies. 

Finally we applied as fiducial limits the p constraints relative to 'MED' and the marginal- 
ized astrophysical background to the neutralino LSP in a simplified version of the pMSSM, 
where we set all the sfermion masses (apart from that of the third generation) to the TeV 
scale. We found that the fiducial Pamela anti-proton and Fermi-LAT gamma-ray limits 
rule out small but non negligible neutralino-chargino mass splittings. In particular for 
m x o < 150 GeV, one can rule out mass splittings up to 20 GeV. Our results also suggest 
that pure wino or wino-like neutralinos are excluded if they are lighter than 450 GeV. 
Overall, this limit surpasses the bounds that can be set by using the XENON100 data and 
even in fact than the projected XENON1T limit. 

Hence from this work, we conclude that present indirect detection data already enable 
one to exclude regions of the parameter space where the neutralino-chargino mass splitting 
is small but non negligible. Since these regions are difficult to probe directly at the LHC, 
these findings show that Fermi-LAT and Pamela data constitute modern tools to explore 
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the supersymmetric parameter space and even beat LHC (and also in fact Direct Detec- 
tion) searches on their own territory, even though - on the negative side - they assume a 
regeneration of the relic density for neutralinos with a very large annihilation cross section. 
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